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Abstract 

The free energy bounds for inhomogeneous current states in doped antiferromagnetic insu- 
lators and spatial configurations of spin and charge degrees of freedom in lightly and heavily 
underdoped phases of this strongly correlated electron system are considered. It is shown that 
states are characterized by a small parameter of geometric origin, which determines the degree 
of packing in the knots of filament manyfolds of the order parameter distributions. We found 
a hierarchy of the spatial scales for flat knots, which gives rise to a free energy gain. A toroid 
moment in an underdoped state is found for the first time. It is shown that in the percola- 
tion picture of charge density distributions, that surround antiferromagnetic nanoclusters, the 
optimal level of the hole density is in qualitative agreement with the experimentally observed 
value. 

PACS numbers: 74.20. Mn, 74.81.-g, 74.78.Na 

1. Introduction 

Doping of planar antiferromagnetic insulators breaks the antiferromagnetic order and gives rise to 
the state that precedes the high-temperature superconducting phase. After one of ten electrons is 
moved, roughly speaking, to the reservoir of dopants, the phase with short-wave length correlations 
of the electron degrees of freedom arises. This strongly correlated electron system has several key 
properties. Among them are strikingly small values of hole doping in the underdoped region of the 
phase diagram, where there are both unusual so-called pseudogap and superconducting phases. 

The "normal" state is characterized by a nonuniform distribution of strongly correlated spin 
P3E] and charge (B1E1Q3E10IE1 degrees of freedom, as well as by an electron-hole asymmetry |9llTfl]. 
A comparison of experimental data obtained from diffent types of superconductors shows that the 
universal two-component dynamics takes place in all cuprates investigated so far. The universal 
multi-component response is clear indication of the presence of the phase stratification pi] in 
cuprates. Study of cooperative behavior in inhomogeneous systems is recognized nowadays as one 
of the most challenging problems of the condensed matter physics |12l \V.\\ and non-Abelian field 
theory ^he aim of this paper is to find bounds on the free energy of inhomogeneous 

current states of doped antiferromagnetic insulators at a moderate level of doping. 

Many recent experiments provided significant support to the idea [TJl that pseudogap phase is 
the initial state for the superconducting phase due to precursor pairing and that it has basically 
the same origin as superconductivity. In this paper we focus on the universal description of both 
phases, considering them on equal footing. Starting from pseudogap phase, we use an equvalent 
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form of the spinor Ginzburg-Landau functional to study a free energy gain. To characterize the 
properties of the spatial configurations of strongly correlated spin and charge degrees of freedom 
in the partial coherent phase we introduce a toroid moment [191 1201 l2~T] . We show how to embed 
considered planar phenomena into three-dimensional space. Using quasi-one-dimensional character 
of spatial distributions of hole density and the percolation picture of the charge configurations, we 
argue also why is the optimal doping in high-temperature superconductors so small. 



2. Model 

To describe spatial distributions of the electron degrees of freedom, we shall use the following 
(3 + 0)1? extended version [23 of the 0(3) nonlinear a model: 

F = F n + F p + F c = (1) 



d 3 x 



\p 2 (d k n) 2 + H 2 k \ + f(d kP ) 2 bp 2 + + f±p 2 c 2 + Fg- 2F lk H tk 



In this equation, the field of the unit vector n(x,y,z) describes spin distributions, while the mo- 
mentum field c = 3/p 2 = a A, where J is the total current, contains paramagnetic (a) and 
diamagnetic (—A) parts. In Eq. (1) Fi k = diC k — d k Ci and Hi k = n • [din x d k n] = did k — d k a,i. For 
simplicity, the Higgs potential in Eq. (1) is expressed only via two constants b and d. The den- 
sity of the current energy yg p 2 c 2 , the surface energy F 2 k and the diamagnetic coupling —2Fi k Hi k 
of the fields a and c determine contributions to the current part, F c , of the free energy. The 
dimensionless units of all the variables are given in Ref. . 

The first term in Eq. (1) is the energy of a Heisenberg antiferromagnet written in the long- wave 
length limit. The self-energy H 2 k of the internal magnetic field Hi = 1/2 Si ks H ks in the long-wave 
length limit [23j includes the mean value < 0|Si[S2 x S3] |0 > of the spin noncollinearity degree 
on the three sites of a plaquette of the square lattice which appears after the electron is moved 
from the fourth site to the reservoir of the dopants. These chiral three-spin correlations of spin 
degrees of freedom are the main building block for the formation of more complicated spin and 
charge structures considered below. The structures are encoded in Eq. (1) by the term (H ik —F ik ) 2 
which describes also mutual frustration of the currents a and c. 

For better understanding of the role of each term in free energy (1), it is useful to compare Eq. 
(1) with the Hamiltonian of the t — J model 

H=^2,Jij SiSj-jmrij - ^j{ a ia a ^ + H.c.) 



<ij>,a 



in the inhomogeneous case. Here ai a = Cj CT (l — n^-a), n. i:a = c^ a c ia , Si = cl a a a pCip/2. This model 
is a standart one when studying different phase states in doped antiferromagnetic insulators. At 
gradient expansion, the first and the second terms in this equation provide the first and the third 
terms in Eq. (1) respectively. The last kinetic term in the t— J model ensures the existence of the 
term (l/16)/9 2 c 2 in Eq. (1). The self-energy B 2 /8tt - curl 2 A - (H ik - F lk ) 2 of the gauge field A 
in Eq. (1) can be also obtained [25} BBj in the framework of the t — J model. 

The exchange integral modulation, the scalar field p(x, y, z) in (1) describes the spin stiffness 
p 2 and depends by the doping level. Indeed, the spin stiffness p 2 relates to the density p\ of 
positively charged vacant sites, so-called "holes", by the constraint p 2 + p\ = 1 [221 - If p 2 1 , 
the hole density p\ , like the term H 2 k , is localized at the boundaries of the regions of the nonzero 
spin density, i.e. it exists on the complement space with respect to the space of the spin density 
definition. Based on this remark and on the lattice patterns of the distributions of the spin degrees 
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of freedom, a low-dimensional analog of the holographic principle, which leads to the conclusion 
that the distribution of the states with broken antiferromagnetic order has the boundary character, 
was proposed in Ref. |28j . 

3. Knots and links of order parameters 

The properties of different phases in model (1) depend on the distributions of competing order 
parameters n, p and c in the ground states. At c = and p = const the configurations of the 

field n(x, y, z) are classified by the Hopf invariant [2J3 1^01 EH Q = j^i J d 3 x Eiki a-i dk a/ . At the 

M 

compactification of M 3 — > M — S 3 and n G S 2 , an integer QeZ^ fl"3(s ) indicates the degree of 
linking or knotting of a filament maniford collection in the three-dimensional space M, where the 
field of the unit vector n(x, y, z) is constant: n(ar, y, z) = n. For two linked loops Q = 1, for the 
trefoil knot Q = 6 and etc. In this case the local minima of the free energy F are determined 32J 
by the inequality: 

F > 32tt 2 |Q| 3/4 . (2) 

The quasi-one-dimensional character of the distribution of the n-field in Eq. (1) at p = const, 
c = 0, if p 7^ const, c = affects also the hole density p\ = 1 — p 2 [2HJ. The configurations of field 
Ph in the form of stripes and loops are most advantageous in this case. 

When the field of the momentum c ^ is unquenched and p = const, free energy bound 
decreases |33j : 

F> 3 27 r 2 |Q| 3 / 4 (l-j^y, (3) 

where the mutual linking number of the fields a and c equals L = J d 3 x Eksi Cfe d s a/ , is 

determined by the invariant L ^ Z. The universal gain of the free energy in Eq. (3) describes the 
condensation energy p2| in the considered model of strongly developed spin distributions. 

The reduction of the free energy bound occurs due to the diamagnetic coupling of the currents in 
Eq. (1). The contribution of this last term is so significant that it consumes completely the energy 
F c and a part of the energy F n . This phase with large pair momentum differs from the FFLO 
state [55, 36 by the unquenched paramagnetic part of the total current. The second distiction 
corresponds in a general case to inhomogenious distribution of the total momentum of pairs existing 
in the incoherent form JHj above the superconducting T c . 

The existence of quasi-one-dimensional configurations with open ends and in the form of loops 
in the framework of the present model, when the free energy density has the form / = (dkp) 2 + 
Veff(p) + (F lk - H lk f, with V eff = -b e ffp 2 + (d/2) p i and b eff = b - (l/16)c 2 - (1/4) (d k n)\ 
corresponds to different phase states [371 138j . In the phase, where the gauge symmetry remains 
unbroken, & e // < and p m in = 0. 

It follows from the equation for the energy V e f / that inhomogeneous distributions of the fields 
exist only in the case b e ff > 0. With the change of doping, approaching the phase with b e ff > 0, 
where the gauge symmetry is broken due to Higgs effect, the fraction of quasi-one-dimensional con- 
figurations with closed ends increases. In the phase b e ff > 0, there exist [38 a quasi-onedimensional 
configurations of length I, satisfying the relation I p m in — 27m. In the ground state c = 0, and 
the value p m in for the amplitude of the Higgs field in the ground state is determined by the con- 
tribution (l/4)(<9fcn) 2 of the characteristic configurations of field n. In this case, distributions in 
the form of a ring are typical configurations of field p j2H] • For greater values (I p m in > 2n) of the 
commensurability parameter I p m in in the ground state, there exist closed quasi-one-dimensional 
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distributions of field p. For sufficiently short (0 < I p m in < 2tt) configurations of the field p in 
the ground state, quasi-one-dimensional field distributions with open ends may be found. Sum- 
marizing, we see that in magneto-inhomogeneous phases, we encounter an internal hierarchical 
multi-scale structure, determined by the density field, which becomes crucial for the description 
of these phases. 

We see also that the state with a finite value of the momentum field c is characterized by 
the competing inhomogeneous order parameters n and c with non-local correlations 1 between 
them due to the entanglement of their one-dimensional distributions. Note that this conclusion 
is in agreement with one of the main statements of Refs. HH] that inhomogeneous current 
states cannot be described with the aid of Ginzburg-Landau functional with the one-component 
order parameter. In our case the existence of inhomogeneous current states is due to the non- 
linear coupling of several order parameters n, c and p at the level of their equations of motion. 
The analogous conclusion may be made, starting from a lattice description of current states in a 
so-called staggered flux phase considering strong correlations of charge and spin degrees of 
freedom in the mean field approach. Flux phases and DDV states i-e. the phases with 

a modulation of the current density, are well-known particular examples of the toroid states. In 
this paper we focus our attention on non-perturbative analysis, and consider a general case when 
all inhomogeneous order parameters are present. 



4. Inhomogeneous current states 

The basic property of model (1) is the compactness of field n, which leads to its one-dimensional 
distributions in the form of knots [22l EDI EH EH ESI S3 > an d to the discrete character of spectrum 
(2). To analyze the general case of c ^ and p ^ const, we first write down the energy F n at 
p = const and x% — > in the form of F n = J d 3 x gi (c^n) 2 + gi H 2 k = Jgigi (r/R + R/r), 

where r is the scale factor, R — \fg2Jg1 and 171 ~ p 2 and gi — 1 are the coupling constants. 
The latter equation shows that the characteristic knot size at the minimum of the energy F n is 
proportional to p . 

With the increase of doping p\, i. e. with decrease of p 2 , the radius R ~ l/p grows up to 
its certain critical value R cr . Beginning from this radius, the configurations of field n, forming 
the knot, are unstable [HI with respect to the scale transformations. This means that knots of a 
size greater than R cr collapse, and form localized structures near the points of the stereographic 
projection R 3 — ► S 3 . This phenomenon appears to be the key one, when we compare different 
contributions to the free energy (1) in an inhomogeneous current state at small enough values of 

We see that at sufficiently small density p 2 or (which is the same) at sufficiently high level of 
doping p\ , beginning from the finite values of density p 2 , in order to describe strongly correlated 
charge and spin degrees of freedom, one may introduce a small parameter, determined by the 
relation Rq/R, where Rq is the correlation length of the n- and c-field distributions. The radius 
i?o has the geometric meaning of the filament thickness, for which the configurations of these fields 
are determined. A typical value of Ro is of the order of 10 A, i. e. it equals three-four lattice 
constants. Since Ginzburg parameter in this case is much less than unity, we may use the mean 

1 In (2 + 1)D dimensional form of the dynamical description, the non-local interaction between electrons may 
be described using Chern-Simons statistical gauge fields a % a , which induce strong correlations of the phases of 
fermion quantum states at large distances. The world line braiding matrix kij in the Chern-Simons action S = 

-P- J d 3 x eap-ya^d/s ai f in this case is proportional to the matrix K ap < (see below), describing the knotting degree 

M 

of the field configurations of the (3 + 0)D free energy. 
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field approach (1). A characteristic scale R has the meaning [12) of the effective radius of current 
correlations. 

Let us assume that Rq < R < R cr and compare the decrease of the free energy (1) due to the 
diamagnetic term —2FikHik with its increase due to the term (dkp) 2 - If the contribution of the 
diamagnetic term —2FikHik is of the order of Cq/Rq, where Co is the amplitude of the momentum 
field c, the estimate of the maximum contribution of the term (dkp) 2 yields 1/(RRq) 2 . The gain 
obviously should be also greater than the contributions of other terms in the last bracket in Eq. 
(1), i. e. of the surface term F? k ~ c \lR\ as well as of the current energy p 2 c 2 ~ F 2 k (R n /R) 2 . 
Combining these inequalities, we see that at a — (Rq/R) 2 << 1 and at a finite value of the 
momentum determined by the inequality 



even in a state with the inhomogeneous distribution of the density p 2 the free energy may be smaller 
due to diamagnetism than at c = in Eq. (1). The parameter a has the meaning of the degree of 
filament packing in the three-dimensional space. It is the searched small parameter in our problem 
of strongly correlated spin and charge degrees of freedom when all coupling constants are of the 
order of unity. In the case n = const, it is proportional to k~ 2 where k is the Ginzburg-Landau 
parameter. 

To find the condition for a state with momentum from the range (4) when the energy is less, 
than at c = (in spite of inhomogeneity) , let us compare the gain of the total free energy with 
the contribution from the term (dkp) 2 . The loss of the energy is of the order of R z /R 2 . Here R z 
is the "vertical" knot size. Therefore, the diamagnetic contribution is the main one for sufficiently 
flat knots [SHI, existing at Rq < R z < R < R cr . This means that the energy of configurations in 
the sector with the topological invariant L ^ even in the inhomogeneous state may be less than 
at c = 0. This conclusion is one of the main results of the paper. The condensation energy 
which depends on the linking degrees L and Q, decreases in comparison with the case of p = const 
y.V,\\ . The minimal value of the free energy gain —AF can be estimated as — AF ~ (647r 2 /<5 1 / 4 ) a. 

4. Toroid state 

Let us assume that the level of doping is small enough, i.e. ph <§C 1. Then the spin density p 2 is 
large and a knot size R is so small that a < 1, The decrease of the degree a of filament packing in 
the knot with the increase of ph is accompanied by the transition to a state, in which among planar 
field projections the field configurations p(x, y) in the form of closed one-dimensional distributions 
are most preferable [28] . 

The fact, whether the boundary charge states at the boundaries of these isulating regions will 
form an infinite percolative cluster, depends on the hole density p\. Let the level of doping be 
insufficient for the appearance of an extended connectable cluster. In accordance with Ginzburg's 
proposal, we will call this phase state with the spontaneous diamagnetism a toroid state [3D|- The 
phase state, in which current order is characterized by a polar vector T, that changes sign under- 
time reversal, can be represented by the ordering of toroid moments. It is characterized by a toroid 
moment T, which is defined in the following way IJH] 



a/Ro < cq < 1/Rq 



(4) 



a = curl curl T . 



(5) 



Using the equality (4) and the identity [H] 



(curia). = - 



e m n[d k n x <9;n] 



(6) 
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one can express the toroid moment via the chiral degree of noncollinearity n • [9tnx din] in the 
form: 



(r-r') fc f 3 „ n • [d k n x djn] 
a r ; — : m 



|r-r'| 3 J I'' - i" 



(7) 



We see that toroid moment T is given by the non-local Biot-Savart law as well as by the factor 
which is determined by the Coulomb Green function. The vector T characterizes the distribution 
of the poloidal component of the current on the torus jlfjl |20] . The magnetic field of this current 
is confined to the interior of the torus. For the degree of noncollinearity (localized at a point) 
n • [9fen x din) = e k i 5(r) of the spin distributions, the toroid moment has the form 



w = S <*v^4> = ^. (8) 



16n 2 J T r ' | r _ r '| 3 8irr 

If the chiral density n • [<9 fc n x c^n] = (e fci /(27ri?)) 8(z)S(r — R) of spin distributions is localized 
in the plane z = on a ring of radius R, the toroidal moment equals li(r, 2 = 0) = Eki{xi/r)t(r) 
with the function t(r) having its maximum at r = R and t(r) ~ r~ 2 at r ^ R. 

In the toroid state the diamagnetic susceptibility x' ls determined ^Oj by the radius R of 
current correlations in such a way that x' — Xl/cy. Here xl = —e 2 kF/(X^ 2 mc 2 ) is the Landau 
diamagnetic susceptibility for non-interacting electrons. For porous knots with a <ti 1, it can be 
rather significant. Anomalously great diamagnetic susceptibility is caused by a great radius R of 
current correlations. It should also be noted that the condition Xl ~ —e 2 nr 2 /(mc 2 ) > —l/Air 
equivalent to the inequality r < ru = c/u> p once again attract our attention to the boundary 
character of the considered phenomena. Here n is the mean electron density, ro is the Debye 
radius, and uj p = (Anne 2 /m) 1 / 2 is the plasma frequency. We pay attention to the fact that in this 
phase the total magnetic moment of clusters equals zero and the toroid state is d- symmetric. 

The realization of the toroid state is achieved in systems belonging to 31 classes of magnetic 
symmetry where the vector T can exist [42]. The conditions for the realization of a toroid 
state are easier to satisfy (i) in the neihborhood of inhomogeneities and (ii) in systems with strong 
electron correlations, (hi) The toroidal state formed against the background of the antiferromag- 
netic state is the most likely situation as well as (iv) conditions of the presence of charged and 
magnetic impurities, i.e. holes and dynamical non-collinear spin distributions in our case, (v) 
Strong spin-orbit interaction also gives rise to effects associated with toroidal ordering. All these 
conditions are met in the system at hand. 



5. Planar phenomena in 3D space 

Under the conditions when all the coupling constants in (1) are of the order of unity 2 the main ob- 
ject of the study appears to be the geometric properties and the topology of knotted one-dimensional 
manifolds of the order parameter domain, which give the main contribution to the free energy. One 
of the results of this approach is the appearance of the parameter a = (Rq/R) 2 < 1, characterizing 
the degree of filament packing in the knot. The existence of the minimal value of the field p 2 , 
playing the role of the running coupling constant in Eq. (1), i. e. the existence of the minimal 
value a m in = (Ro/R cr ) 2 , is similar to the emergence of the solution for the superconducting gap in 
a two-component model [52j, beginning (unlike BCS model) with the finite value of the coupling 
constant. 

Let us describe a method of embedding planar phenomena into the 3D space. We choose a 
tuning parameter p 2 ~ 1/i? 2 that a m i„ < a << 1. Let us consider the projection of complementary 

2 This is one of the reasons, why the problem discussed cannot be currently solved, even numerically. 
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Figure 1: (colour online). Complementary surfaces for a flat (l.h.s.) and extended along z-axis 
(r.h.s.) trefoil knots. 

spaces for knots which are pulled over the knots (so-called Seifert surfaces) to the plane (x, y). For 
definiteness, we consider the simplest (among torus knots) trefoil knot with Q = 6. In the case of 
the flat (R z < R) knot, the complementary space is a triple twisted band. In the case of the knot 
stretched along the axis z (at R < R z < R cr ) the Seifert surface is a torus with a hole (see Fig. 1) 

3 

The charge degrees of freedom are distributed (as it has been already mentioned) at the bound- 
ary of the complementary space for knots. Therefore, the boundary line presents a circle, which 
is the projection of these one-dimensional boundary configurations for a sufficiently flat knot. For 
the stretched knot (at R < R z < R cr ), the line is a stripe with the length of the order of Rq, equal 
to the scale of the hole mentioned above. Summarizing, we can see that the free energy gain in an 
inhomogeneous current state is achieved, due to diamagnetism, for the distributions of the charge 
density on the plane in the form of loops. 

The study of such configurations of the charge density is based on measuring the density of 
states (DOS) by means of scanning tunneling microscope (STM). In such experiments jSH] the 
formation of the DOS filaments in the form of rings have been discovered recently in the entan- 
gled labyrinthine pattern of quasi-one-dimensional distributions of the charge density in planar 
compounds Ca2- x Na x CuC>2Cl2 at the moderate doping x. 

6. Optimal doping 

Whether the phase both with a cr < a << 1 and the finite pair momentum l-'ifij Co remains su- 
perconducting or appears to be an intermediate state with spontaneous diamagnetism 0^1 depends 
perhaps on the connectness of the total current cluster (see e.g. [l4|), formed by the projections 
of Seifert surfaces on the plane. 

It is natural to suppose that the occurrence of pre-formed electron pairs in the underdoped state 
[18] and a phase coherence in the superconducting phase are indepent phenomena. Starting with 
the picture of pre-formed pairs, we discuss in this section some features of the phase coherence 
due to percolation |SS] using charge density structures mentioned above. With the increase of 

3 The braided filaments of the trefoil, stretched along the time axis, one can consider in (2 + 1)D dynamical case 
as the result of the action of bf on the quantum states of the braid group operator b\ . 
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the hole density p\, the size of ring-shaped current distributions increases to such an extent that 
the conditions arise for the appearance of two-dimensional connectable current clusters and the 
transition to the state with c ^ with the gain of free energy. We shall find the level Pbest, for 
which the temperature of the transition to the superconducting state is maximal. 

If we did not take into account the distributions of the charge degrees of freedom at the 
boundaries of the filament configurations of the field n in the knot only, we would deal with 2D 
percolation problem with a well-known p2j answer p cr = 0.5 for the bond percolation threshold 
value. A considerable distinction between the experimentally observed value Pbest, ex P = 0.16 from 
this value is a convenient test for checking the accepted point of view. 

To solve the problem in question we shall use the recent experimental results - With 
the aid of the scanning tunnel microscope it was found in j^Hj that among quasi-one-dimensional 
configurations with an increased density of charge states in the underdoped region, the fraction 
of filament field distributions in the form of ring-shaped structures increases with the increase 
of the doping level. Using the corresponding trial functions of paper j2H] it was shown that in 
the framework of model (1), such configurations of the field pt are preferable when the degree 
of doping increases. The interpretation of the results of measuring the distributions of the spin 
degrees of freedom in underdoped states on the basis of neutron scattering p][2] is also based on the 
assumption of quasi-one-dimensional nanoscale distributions in the planes of conducting clusters 
in the form of thin rings. 

Let us suppose that the boundary currents are distributed along the boundaries of ring insulator 
regions with such a relationship that the ring thickness differs from its diameter, as it follows from 
the experiment by one order (a 1 / 2 ~ 0.1). The geometrical confinement allows us to neglect 
the contribution of internal (with respect to ring-shaped) distributions of degrees of freedom. Thus 
for the threshold of the flow along the connectable cluster, formed by ideal distributions in the 
form of rings, we get a considerabl smaller value: Pbest — 0.5 — 0.314 = 0.186 in comparison with 
the standart one p cr = 0.5. For a current cell in the form of a hexahedron, Pbest = 0.179. This 
example shows that comparatively slight deviations of the form of the elementary current cell 
from the ideal ring-shaped form may lead to a shift in the value Pbest towards the experimentally 
observed value Pbest, exp = 0.16. If a square, constructed by the distributions in the form of stripes, 
were used instead of a ring, we would have pbest = 0.5 — 0.4 = 0.1. 

Study of general properties of high-temperature superconducting states, belonging to different 
classes of compaunds, shows that the dependence of the critical temperature on the hole concen- 
tration, for example in Y\- x Ca x BaiCuzO§ and La<i- x Sr x Cu)±, is well represented [561 I57j by 
the empiric curve T c /T CtJnax = 1 — 82.6(p — 0.16) 2 at p — x/2, This parabolic dependence T c (p), 
has zeros at p m in ~ 0.05, p ma x ~ 0.27 and a maximum at pbest — 0.16. Optimally doped com- 
pounds based on thallium, for example TZ0.5P&0. 5^2^00. 8^0. 2C , M2C ) 7, have Pbest = 0.15. Some 
other classes of high-temperature superconductors point to similar values. We note also that the 
underdoped pseudo-gap phase with strong developed spin and charge distributions appears in the 
above mentioned classes of the compounds at a well-defined p>H] hole concentration p w 0.19. 



6. Conclusion 

Under conditions of moderate doping p\, the value of p 2 may be such that the knot size R ~ 1/p ~ 
R z appears to be comparable with the correlation length Rq. Since all values in Eq. (1) are of 
the same order, these greater values of p 2 correspond to smaller values of the momentum c and 
due to Eq. (3) the energy gain is less. In the optimum case of great values c the self-dual state 
F n ~ F c , which is characterized by the dense packing of knot filaments (a < 1), a— c = A = and 
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by the matrix jS3 K a/3 = J d 3 x e lk iafd k af = y L Q' ) ^ ® \ 1 1 J' where a i = 

af = Ci, may be considered as initial compressible state for the transition to the incompressible 
phase a << 1 with the increase of the doping degree. The compressibility of the state at a < 1 
occurs due to the existence of non-zero Poisson bracket {x ,y } for the coordinates x ,y of the 
vortex filament center. The uncertainty of the location of the core center at the dense filament 
packing gives rise to the nullification of the linking degree Q due to the effect of the filament 
reconnection. This leads to the transition to the ground state with Q = 0. 

The comparison of the results of the present paper with the conclusions of |5Qj, following from 
the application of the conformal field theory methods to find the exact solution of the problem 
of superconducting states in cluster systems, requires thorough analysis of the details of the nul- 
lification (fc + 2 = 0) mechanism of the linking coefficients k + 2 ~ Q in SU (2) Chern-Simons 
action. 

In conclusion, we have compared the value of the free energy gain due to diamagnetism in the 
inhomogeneous current state with the value of the energy loss due to inhomogeneity and found the 
conditions, under which the energy of inhomogeneous states with a current appears to be less than 
in the absence of the current. Studying the threshold level of doping with the use of the percolation 
approach of describing current configurations, existing in doped antiferromagnetic insulators, we 
have found, that in the toroid state the optimal value of the hole concentration is close to the 
experimentally observed value. The toroid moment in this d-symmetric state is expressed via the 
chiral degree of noncollinearity of the distribution of the spin degrees of freedom. 
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